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Abstract 

We show that a signal can propagate in a particular direction 
through a model random medium regardless of the precise state of 
the medium. As a prototype, we consider a point particle moving on 
a one-dimensional lattice whose sites are occupied by scatterers with 
the following properties: (i) the state of each site is defined by its spin 
(up or down); (ii) the particle arriving at a site is scattered forward 
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(backward) if the spin is up (down); (iii) the state of the site is mod- 
ified by the passage of the particle, i.e. the spin of the site where a 
scattering has taken place, flips (T^i)- We consider one dimensional 
and triangular lattices, for which we give a microscopic description of 
the dynamics, prove the propagation of a particle through the scatter- 
ers, and compute analytically its statistical properties. In particular 
we prove that, in one dimension, the average propagation velocity is 
(c(9)) = 1/(3 — 2q), with q the probability that a site has a spin t, 
and, in the triangular lattice, the average propagation velocity is in- 
dependent of the scatterers distribution: (c) = 1/8. In both cases, 
the origin of the propagation is a blocking mechanism, restricting the 
motion of the particle in the direction opposite to the ultimate prop- 
agation direction, and there is a specific re-organization of the spins 
after the passage of the particle. A detailed mathematical analysis of 
this phenomenon is, to the best of our knowledge, presented here for 
the first time. 



1 Motivation 

Discrete systems with simple microscopic dynamics can exhibit peculiar be- 
havior showing some degree of complexity on a large scale. These systems 
have raised particular interest because they can be viewed as paradigms for 
complex phenomena such as growth processes, signal propagation in ran- 
dom media, spatio-temporal structuring in excitable media, or evolution- 
ary dynamics In particular, phenomena such as anomalous diffusion 
and oscillatory propagation have been reported by Cohen and co-workers 
[0] who studied, numerically and theoretically, the motion of a particle in 
two-dimensional Lorentz Lattice Gases whose sites are occupied by scatter- 
ers which deflect the particle according to an a priori given rule. The most 
striking behavior in the particle dynamics is observed when the rule includes 
a feed-back of the particle on the substrate: the passage of the particle mod- 
ifies the scattering property of the visited sites, e.g. if a particle arriving 
at a site is scattered say to the right (R), the state of the site is changed 
(R =^ L) such that on its next visit the particle will be scattered to the 
left ("flipping scatterers"). This "interaction" between the particle and the 
state of the scatterer modifies considerably the dynamics as compared to the 
classical Lorentz Lattice Gas. For instance, lattices with flipping scatterers 



2 



can yield oscillatory motion with overall propagation [Q, while random fixed 
scatterers produce closed trajectories. Bunimovich and Troubetzkoy proved 
a number of theorems about the boundedness or unboundedness of the tra- 
jectories of the particle for certain lattice models, and also discussed some 
generalizations 0. 

In order to illustrate the propagative behavior, we show in Fig.l the 
case of a random Delaunay lattice 0: here the particle arriving at a node 
is defiected with the largest possible angle, either to the right or to the 
left, depending on the state of the scatterer (which after the passage of the 
particle goes into the reverse state); whatever the initial configuration of the 
scatterers on the lattice, the particle will, after a few time steps, always enter 
a propagation phase (see the propagation strip shown in Fig.l). 

Here we present a detailed mathematical analysis of unbounded trajecto- 
ries in one-dimensional and two-dimensional (triangular) lattices. In partic- 
ular we show that for any initial configuration of the scatterers, the particle 
will always propagate in a certain direction with a given average velocity. 
The direction of the propagation is determined by the initial velocity of the 
particle and the initial state of the scatterers in a small neighborhood of the 
initial position of the particle. On the triangular lattice the entire trajec- 
tory becomes quickly confined to a particular strip, which is bounded by two 
adjacent parallel lines of the lattice 0. The propagation is due to a blocking 
mechanism which prevents the particle from moving in a direction opposite 
to the propagation direction for more than a few steps; thus the particle can 
visit any lattice site at most three times in one dimension, and six times on 
the triangular lattice. The blocking mechanism is due to an organization of 
the scatterers induced by their interaction with the moving particle, which it- 
self is thereby forced to propagate indefinitely in a particular direction. While 
in one dimension the rearrangement of the scatterers after the passage of the 
particle is a simple mapping of the initial configuration, on the triangular 
lattice a self-organization of the scatterers takes place along the propagation 
strip: the visited sites on one of the two boundary lines are occupied by only 
one kind (R or L) of scatterers, while the other boundary line is occupied 
by only the other kind (L or R). We find another important distinction: the 
average propagation velocity (averaged over all possible initial distributions 

^This propagation phenomenon was first observed and analyzed on the triangular lattice 
by Kong and Cohen 
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of the scatterers with a fixed ratio of R and L scatterers) in one dimension, 
depends on the a priori probabihty that a site has a spin |, while in the 
two-dimensional (triangular) lattice, it is a constant (independently of the 
probability that a site has a L or R scatterer). 

We first consider the one dimensional lattice. In section H, we define 
the 1-D system considered here, introduce the basic equations describing 
the microscopic dynamics, and prove propagation of the particle; section ^is 
devoted to a detailed study of the statistical properties of the particle motion. 
The triangular lattice is treated in section ^ where particle propagation and 
its general and statistical properties are analyzed. All our analytical results 
are shown to be in agreement with our corresponding numerical simulation 
data obtained for one- and two-dimensional lattices. We conclude with some 
general comments (section |^) . 

2 The 1-D lattice 

We consider a regular one- dimensional lattice where a single particle moves 
from one site to the next, with the following properties: (i) the state of each 
site is defined by its spin (up or down); (ii) the particle arriving at a site is 
scattered forward (backward) if the spin is up (down); (iii) the spin of the site 
where a scattering has taken place, fiips (T^i, i^T) Q- Setting the distance 
between neighboring sites and the speed of the particle equal to unity, the 
particle moves from one lattice site to the next in one unit time step . The 
initial position and velocity of the particle are arbitrarily fixed, and the spins 
are arbitrarily distributed on the lattice with an a priori probability q that 
a site has a spin ]. 

2.1 Definitions 

(i) Define the position R{t) of the particle at time t as the location of the 
site r where the particle resides at time t. 

generalization to more complicated transitions in the scattering properties of the 
sites has been considered for cyclic cellular automata by Bunimovich and Troubetzkoy ^ 
who proved a number of general theorems, without having to use explicitly the equations 
of motion for the particle and the equations for the dynamics of the scatterers. 
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(ii) Define the velocity of the particle: C{t), whose value is + 1 (when the 
particle moves to the right) or — 1 (when it moves to the left) 0. 

(iii) Define the spin (the orientation of the scatterer) Q at site r, at time t, as 
a Boolean variable: r]{r,t), with value + 1 {spin up, which does not modify 
the velocity of the particle), or — 1 {spin down , which reverses the velocity 
of the particle from C = + ltoC = — lor vice versa). 

2.2 Basic equations 

(1) At each time step the particle moves to a neighboring site according to 
its velocity, and its new position is given by the equation of motion: 

R{t + 1) = R{t) + C{t). (2.1) 

(2) The velocity of the particle at its new position is maintained or reversed, 
depending on the state of the spin (the orientation of the scatterer): 

C{t + 1) = C{t)r]{R{t + l),t). (2.2) 

In ri{R{t + l),t), t refers to the time corresponding to the state of the spin 
before flipping; we shall denote by t+ the time just after the particle arrived 
at a site and made its spin flip (t^^i)- 

(3) The spin of the site hit by the particle is reversed (all other spins remain 
unchanged) : 

7]{r, t + 1) = r]{r, t) (1 - 5r,R(t+i)) - vi^, t) 6r,R{t+i) 

= v{r,t){l - 2Sr,Rit+i)); (2.3) 

so 

r]{R{t + 1), t + 1) = - r]{R{t + 1), t) . (2.4) 

Equations {\i-l\), (^^.2]), and ( f2.4D are the basic equations of the dynamics 
on the 1-D lattice. 

■^With the convention that the positive direction on the 1-D lattice is from left to right, 
we can omit vectorial notation for simplicity. 

^In previous publications the scatterers were indicated as to their right or left scat- 
tering properties (R or L); for the 1-D case a representation of the scatterers by up or 
down spins is more convenient. 
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2.3 Propagative motion 

We say that the particle propagates in one direction on the one-dimensional 
lattice, if it visits any lattice site not more than some fixed finite number 
of times. This definition implies that the trajectory of the particle tends 
to +00 or —00, when the time goes to infinity. Indeed since the particle 
cannot return to a site after it has visited it a finite number of times, the 
particle must remain on a semi-line to the right (or to the left) of this site, 
which holds equally for the next visited site. Consequently the position of 
the particle must approach -\-oo or —00 when t ^ 00. 
We now formulate the propagation theorem. 

Theorem 1 : A particle moving from site to site in a one-dimensional lattice 
fully occupied with fiipping scatterers (spins), propagates in one direction, 
independently of the initial distribution of the spins on the lattice. The prop- 
agation direction depends upon the orientation of the spin at the origin of 
the particle motion (r = 0) and of that at the site r — +1, if the initial 
particle velocity C(t = 0+) = C{0+) — +1, or alternatively of that at site 
r = -1, if C(0+) = -1. 

Proof : Without loss of generahty we assume that the initial velocity of the 
particle is positive, i.e., C{0+) — 1. Then we have to consider two cases 
according to whether the spin of the site at the origin {r — 0,t — 0) is either 
t or I (see Fig. 2 for illustration of the proof). 

A. If at t = 0+, ri{R{0), 0+) = r]{0, 0+) = —1, there are two possibilities: 

(1) the spin at time at the position R{1) is t , i.e. ri{R{l), 0) = r]{l, 0) = 1. 
Then at time 1+, C(l+) = 1 and rj{R{l), 1+) = rj{l, 1+) = -1; 

(2) the spin at time at the position R{1) is | , i.e., ?7(-R(l),0) = 77(1,0) = 
-1. At time 1+, C(l+) = -1 and 'n{l,l+) = +1, while at time 2+ , 
C(2+) = +1 and r/(i?(2),2+) = r^(0,2+) = 1. Then at time 3+, C{3+) = 1 
and r;(i?(3), 3-I-) = r/(l, 3+) = -1. 

Therefore in both cases (1) and (2), the system is in a state where the 
particle is at r = 1, leaving the visited sites with the velocity +1 (C(l+) = 1 
in case (1) and C(3+) = 1 in case (2)), while the scatterer at this site (r = 1) 
is I . Thus we retrieve the initial situation shifted by one lattice unit to the 
right, and consequently the same analysis can be repeated for the next site, 
and so on. 
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We call the state of the system when C{t+) = 1 and r]{R(t),t+) = —1, a 
left blocking pattern. So the above analysis shows that a left blocking pattern 
at the site R{t) will be shifted to a left blocking pattern at the next site to 
the right oi R{t) at R{t) + 1, in one time step in case (1) + = R{t) + 1), 
or in 3 time steps in case (2) {R{t + 3) = R{t) + 1), leading, in both cases, 
to propagation of the particle to the right. 

B. If at t = 0+, r]{R{0),0+) = ?7(0,0+) = +1, there are again two 
possibilities: 

(1) the spin at time at the position R{1) is | , i.e. ri{R{l), 0) = ri{l, 0) = +1. 
Then C(l+) = 1 and ri{R{l), 1+) = ^^(l, 1+) = —1, and the same blocking 
pattern obtains at r = 1 at the time t = 1+, as in case A(l). 

(2) the spin at time at the position -R(l) is J, , i.e. ri{l, 0) = —1. Therefore 
C(l+) = -1 and r7(i?(l), 1+) = 7^(1, 1+) = 1. Then at time 2+, C(2+) = -1 
and rj{R(2), 2+) = 77(0, 2+) = —1, and we obtain a situation opposite to that 
in Case A(2), which implies that there is now a n^f/ii blocking pattern leading 
to a propagation to the left. Thus the same analysis is applicable here, if 
one changes positive velocities and position changes to negative ones and 
vice- versa. 

We conclude (see Fig. 2) that there are three cases where the initial con- 
ditions for C and 77 produce a left blocking mechanism and a propagation to 
the right: 

A(l): C(0+) = +1, 77(0, 0+) = -1, 77(1, 0) = +1, 
A(2): C(0+) = +1, 77(0,0+) = -1, 77(1,0) = -1, 
B(l): C(0+) = +1, 77(0, 0+) = +1, 7^(1, 0) = +1, 

and one initial condition: 

B(2): C(0+) = +1, 77(0,0+) = +1, 77(1,0) = -1, 

which produces a right blocking mechanism and a propagation to the left . 

Clearly the same analysis holds if the initial velocity of the particle is to 
the left: for C(0+) = — 1, there are three cases where the propagation is to 
the left, and one where the propagation is to the right. Propagation in the 
direction of the initial velocity occurs in three out of four cases. Notice that 
by fixing the initial condition with only one spin: r]{R{0), 0+) = —1, i.e. the 
spin is set J, at the origin, the particle always propagates in the direction of 
its initial velocity. 
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Two corollaries follow from Theorem 1. 
Corollary 1 : A particle moving on a 1-D lattice fully occupied with flipping 
scatterers visits any lattice site at most three times. 

The proof follows immediately from Theorem 1 asa consequence of the block- 
ing mechanism. 

Corollary 2 : The propagation velocity between two sites separated by a 
distance i containing u spins | , is Cg = {3 — 2u/i)^^. 

Proof : We first note that the time taken by the particle to cover the 
distance i equals i plus twice the number of times the particle visits a site 
with spin J, , forcing the particle into a backward-forward motion. Thus 
ti = i + 2{i — u) = 3i — 2u. Since Ci = i/ti, the corollary follows immedi- 
ately. 

Finally we note that after the passage of the particle through the lat- 
tice, all spins of the visited sites have been reversed (t"^i) with a shift of 
one lattice site opposite to the ultimate velocity of propagation. This re- 
organization of the spins, which is a consequence of the blocking mechanism, 
is illustrated in Fig. 3. 

3 Statistical properties 

Statistical properties are obtained by considering average quantities (•••), 
where the average is taken over all initial spin configurations, with the initial 
a priori probability q that a site is in the state t- 

3.1 Average velocity 

From Corollary 2, it follows immediately that the average time (te) taken by 
the particle to cover a finite segment of length £ on the 1-D lattice is: (ti) = 
3i — 2{u) = £{3 — 2q). Then the average propagation velocity (c(g)) = c{q) 
is given by 

c{q)^e/{te)^l/{3-2q), (3.1) 

independently of the direction of propagation and regardless of the way the 
spins are organized. Examples are shown in Fig.4. 
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3.2 Time evolution 



The equation of motion of the particle can be written in terms of two Boolean 
variables: 

(i) the occupation variable n{r, t), which is equal to 1 if the particle is at site 
r at time t for the first time; otherwise n{r,t) = 0; 

(ii) the state of the spin at site r at time t: ^^(r, t) = 1 if the spin is t , and 
^j(r, t) = 1 if the spin is | , with Sj=|^j = 1.0 Then the equation for n(r, t) 
reads 

n(r, t) = ^^(r - 1, 0) n(r - 1, t - 1) + ^^{r - 1, 0) n(r - 1, t - 3) , (3.2) 

which is a microscopic equation. 

We now define Pi(r, t) = {n{r, t)), the probability that the particle, start- 
ing at the origin, visits site r at time t for the first time. In the following 
we shall use f{r,t) = Pi(r, t) for short, with the property I]t-Pi(r, t) = 
Y^f- f{r, t) = 1, V r > 0, which follows from the fact that each site will have 
a first visit. We also have (^^(r, 0)) = q and (^j,(r, 0)) = 1 — g, V r. Averaging 
the basic equation ( |3.2| ) over all initial configurations and using the fact that 
{^n) on the r.h.s. of (|3.2|) can be written as {C,){n) then yields 

/(r + l,t + l) = qf{r,t) + (1 - g) /(r, t - 2) , (3.3) 

which shows a formal analogy with the equation for the biased random walk 
(BRW)0. For convenience in the forthcoming development, we perform a shift 
of variable (t — > t + 3) to rewrite Eq.(p.3|) as 



/(r + l,t + 3) = g/(r,t + 2) + (l-g)/(r,t). (3.4) 

Considering that the particle, arriving from the left, is at r = for the first 
time at t = with unknown velocity, and that the spin at r = — 1 is down, 
the initial conditions are given by 

/(r,t = 0) = 6r,o; /(r,t = l) = /(r,t = 2) = 6r,+2 ■ (3.5) 



^Notice the difference between rj = ±1 in section |2.2| ||] and = {0, 1}, which is the 
quantity used in practice for numerical simulations. 

^The BRW equation § reads g(r+ 1, t+ 1) = hg{r,t) + {l~h) g{r + 2,t) , h ^0,1, 
where g{r,t) is the distribution of the random walker which, at each time step, moves to 
the right or to the left with probability h and 1 — h respectively (for h = 1/2, the equation 
describes the usual random walk). For a comparison with Eq.(3.3), see also Q]. 
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Introducing the space-Fourier transform fk{t) = J^{f{r,t)} and its discrete 
Laplace transform fk{s) = C{fk(t)}, we note that 

CJ^{f{r + n,t + m)} = e'>'^[e"'^ fk{s) - e"'' MO) - e^"'-'^ Ml) - 

eJkim-1)]. (3.6) 



Equation ( p.4| ) is then Fourier-Laplace transformed to yield 

e^V [e'[e'{fk{s) - MO)) - /,(!)] - M'2)_ = 

qe^[e^{Ms)-MO))-M^)] + {1 - q) Ms) , (3.7) 

with 

/,(0) = 1; /,(!) = ge'^ M'^) = e'^' , (3.8) 

which follow from the initial conditions (|3.5|) . Setting z = e'^ and k = e*'^, 
the solution to Eq.(3.6) reads 

f(,-. _ z^ + q{K~l/K) + q\K'-l) 

' ~ z\z - qlK) - (1 - q)/K • ^^-^^ 

The function fniz) can be explicitly inverted in time and space for q = 1 
and g = 0; the results for f{r,t) = Pi{r,t) are given, as expected, in terms 
of (5-functions, S{r,t) for g = 1 and S{r,t/3) for g = 0. When g 7^ or 1, we 
consider the case k = 1, i.e. the limit of long wavelengths [k = 0), for which 
the poles of the above transform can be computed easily []; we obtain 



L=i{z) = z^ [{z - zo){z - z+){z - 2_)]-i , (3.10) 



where 



with 



z± = (1 — g)"^/^(costi7 =1= « sintxj) , (3.11) 



cosw = --(1 -g)^/^ sinro = i (3 + g)^/^ (3.12) 



^The computation for /t ^ 1 (with q 7^ or 1) is quite involved and the exphcit 
analytical results are not very instructive. 
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The function /k=i(-2) is then inverted (using residues) to obtain the exphcit 
expression for fk=o(t) = T,rf{r,t); after some straightforward algebra, we 
find 

fk=o{t) = [1 + 2(1 - q) e-«* (cos wt + ^ sin wt)] , (3.13) 

with 



C = - log , = arctan — = -, zri^, ttttt- (3.14) 

^ 2^1-g' VI- ? (1 -g)i/2(3 + g)V2 ^ ' 

The case g = provides a trivial but illustrative example ; it is straight- 
forwardly verified that the solution then reads 

fl~,{t) = ^(l + 2cosyt), (3.15) 

which is a repeated sequence of values (1,0,0) as can be inferred from the 
dynamics of the particle when all spins are initially down. The case q = 1 
(all spins initially up) yields the obvious result /|ro(t) = 1. 

The long-time behavior of fk=o{t) is easily obtained from the limit s — > 
of fk^Q^s) by expanding fk=o{s) to lowest order in s 

lim fu-n(s) = , (3.16) 

which yields by Laplace inversion 

hm fk=o{t) = . (3.17) 

S — 2q 



This result, combined with ( p.lj) , shows that 

lim ^Pi(r,t) = c(g). (3.18) 



t— >oo ' 

r 



In Fig. 5 we show that the above analytical results and the simulation data 
are in perfect agreement. An alternative derivation of ( |3.18|) , directly based 
on the analysis of the particle dynamics, is given in Appendix A. 
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3.3 Space evolution 

We now construct an explicit analytical expression in r and t for the prob- 
ability that the particle visits, for the first time, site r at time t, i.e. the 
function Pi{r,t). Consider, on the 1-D lattice, a segment with length i (in 
lattice units) and with spin configuration {r/j} = {rjo, rji, - ■ ■ , rji), where r/j = t 
or [. Suppose that at time zero, the particle is at site r = with velocity 
C(0+) = +1, and the spin is t]q = |. Then the time r taken by the particle 
to travel the distance £ is T{i, {rji}) = 1 + iV| + 3iVj, where N-^ and Ni are the 
numbers of sites with spin up and spin down, respectively, between r = 1 to 
r = i — 1 (since the spin rji is unimportant). So = i — 1 — N-^ and 

T{i,{7]}) =3i-2{l + N^). (3.19) 

The probability for the particle to be for the first time at site r at time t for 
a given spin configuration is 

Pi{r, t; {r]i}) = 5t,r(r,{,,,}), (3.20) 

so that 

P,{r,t) = J2St,rir,{,^})Pi{Vi}), (3.21) 
iVi} 

where P{{r]i}) is the probability of the spin configuration {rji} 

P{{V^})=( "~^]q''^{l-qr-'-''^ , (3.22) 



with q the probability that a site has spin |. Furthermore, using ( |3.19|) for 

i = r, we have 



r-1 / -I 

r — 1 



N^=o \ T / 



|(3r-t-2) 



r — 1 
Ut-r) 



(3.23) 



2 



12 



and it is a matter of simple algebra to show that ( p.23| ) is a solution of the 
difference equation (|3.3| ). This solution is valid for any spin configuration 
with A^i = 0, ■ ■ ■ , r — 1. For the trivial cases where all spins are either up or 
down, one can easily verify the obvious results 

Pi'^"~\r, t = r) = l; Pi'^\r, t = 3r - 2) = 1 . (3.24) 



An example of space-time evolution of the system based on Eq.( |3.23| ) is shown 
in Fig. 6. 



3.4 Long-time large-distance behavior 

An interesting result follows from the computation of the average time the 



particle takes to cover a distance r when r is large. In ( |3.23| ), we set (1— (?) = p 
and {t — r) = 2a, and we consider r large (i.e. r ^ 1); then Pi{r,t) can be 
rewritten as 



Piir,a) 



(3.25) 



Using a standard procedure of probability theory p , we have 



(a) 



E 



q=l-p 



9 , .r 

Pg^ip + q) 



rp,(3.26) 



q=l-p 



which we rewrite as 
1 



it) 



((t) -r) = r(l-g), or 
2 r 



3-2g, 



(3.27) 



i.e., we retrieve the expression for the average velocity c{q) = (3 — 2q)^^. 

We now compute the explicit analytical expression of Pi{r,t) for large r. 
We consider p, q fixed, and r ^ 1, such that -|a— rp| 0, which is equivalent 
to [- — (3 — 2q)] — »• 0, the limit we expect for long times. Then using Stir- 
ling's formula in the binomial coefficient of ( |3.25| ), and performing a Taylor 
expansion in 6a = a — rp, |^ we obtain 



Pi{r,a) 



|r>l 



{2'nrpq) exp(— 5a^/2rpg) . 



(3.28) 
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Noting that 5a /r = ^[(t/r) — (3 — 2g)], ( p.28| ) can be rewritten as 

1 (r - {r{t))r 



with 

(r(t)) = t(3 - 2q)-' 



'27T{rpqy/' 



exp 



2'~fr 



(3.29) 



c{q)t , and 7 = Apq{3 - 2q)-^ = Apqc\q). (3.30) 



In Fig. 7 we show a comparison between simulation data, the computation of 
the binomial expression ( p.23| ), and the analytical result ( p . 29|) . 

The connection with the results of subsection p.2| , in particular with the 
function fk=o{t oo), is easily established by taking the Fourier transform 
of ( |3.29| ) in the limit k = 0, that is dr Pi{r,t)] by setting r = x'^,a = 
(3 — 2g)/y^8pg, and P = a{r) , we obtain from (|3.29| ) 



1 f°° 

- / dxPi{x) 
I Jo 



1 



2V2^(pg)i/2 
which is exactly the result ( |3.17|) .p| 



3-2g 



(3.31) 



3.5 Total probability distribution 

In section |2]^ we showed that any lattice site can be visited at most three 
times {Corollary 1). So we define Pi{r,t), P2{r,t), and P3(r, t), as the proba- 
bilities that the particle be at site r at time t for the first, second, and third 
time, respectively We also introduce the probability P2i{r, t) for the particle 
to be at site r at time t for the second time, without ever having visited the 
next site (at r + 1); then the particle must have visited site r two time-steps 
earlier and this was a first-time visit and the spin was down (see Fig. 8a). 
Therefore 



P2i{r,t) = (l-g)Pi(r,t-2) 



(3.32) 



Equation ( ^.3] ) for /(r, t) = Pi{r, t) is regained by noticing that a first visit at 
site r -|- 1 at time t -|- 1 results from either a direct displacement of the particle 



^Note that the factor 1/2 on the Lh.s. of ( |3.3l| ) comes from the fact that the binomial 
coefficient in (3.25) must be interpreted as zero if r and a are not of the same parity since 
r and t must have same parity. 
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from r to r + 1 when the spin at r is up, i.e. qPi{r, t), or a second visit at r 
two time-steps earher, as described above; so the equation for Pi (r, t) reads 

Pi(r + l,t + l) = qPi{r,t) + P2iir,t) 

= qPi{r,t) + (l-g)Pi(r,t-2), (3.33) 

which is exactly Eq.( p.3| ), with f{r,t) = Pi(r, t). 

Now to obtain an expression for P2{r,t), two cases must be considered: 
the first is that of Fig. 8a, which yields the contribution given by ( ^.32| ), and 
the second, illustrated in Fig. 8b, contributes g(l— g)Pi(r, t—2). Consequently 

P2(r,t) = (l-g2)p,(r,t-2). (3.34) 

There is only one situation which produces three consecutive visits to the 
same site, as shown in Fig.8c: the site must have been visited for the first time 
four time-steps earlier, and the particle must have moved twice backwards 
during these four time-steps, which yields 

P3(r,t) = (l-g)^Pi(r,t-4). (3.35) 

So the total probability that the particle be at site r at time t reads 

P(r,t) = P,(r,t) + il-q')Pi{r,t-2) + il-qfPiir,t-A), (3.36) 

with the property J2r ^) = 1) 

The long-time behavior of the total probability P(r, t) is readily obtained 
by taking t 3> 1 in ( ^.3(j| ), which yields 

P{r,t)\t»i = (3-2g)Pi(r,t)|,»i. (3.37) 



Since t ^ 1 also implies r ^ 1, it follows from ( |3.29| ) that 

where 7 = ic^pq, {r(t)) = c{q)t (see ( p.30D ). Alternatively we write ( |3.38| ) as 
, [2 3-2q it-(t)Y 
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with (t) = r/c{q). Then the following properties, obtained from ( |3.38| ) and 

(EH), 



P(. = --) = ^g^)"'^. (3.40) 



and 

{{t-{t))Y^^ = 2^qr^/^ (3.41) 

show that , for long times, the amplitude of the probability distribution 
P{r,t) decays like and its width grows like ^/r. 

We also note that, by taking the sum over r in ( p.36| ), followed by the 
long-time limit (t 3> 1), we obtain 

1 = (3-2g) lim VPi(r,t), or limVPi(r,t) = c(g) , (3.42) 

t — ^oo — ^ t — ^oo — ^ 

r r 

in accordance with our previous results (see subsection |3.2| ). In Fig. 7 we 
show the long-time behavior of Pi{r, t), which is the same (within a constant 
factor, see ( p. 37] )) as that of P{r,t). 

Finally it is interesting to note that P{r,t), as given by ( |3.39| ), is the 
solution of the differential equation 

drP{r, t) + - dtP{r, t) = I (Apq) d^P{r, t) , (3.43) 

which can be shown to be the continuous limit of the difference equation for 
P{r,t) 0. The appearance of a Gaussian distribution in time, ( p. 391) , and 
the differential equation with propagation ( |3.43D for long times are natural 
consequences of the averaging over the initial spin distributions. 



4 The 2-D triangular lattice 

We consider the motion of a particle on a triangular lattice fully occupied by 
rotators (flipping scatterers) which rotate the velocity vector of the particle 
by an angle of 27r/3 either to the right or to the left, depending on whether 



16 



the rotator on the site hit by the particle is a R or a L scatterer, respec- 
tively.0 The state of the rotator at the site where the scattering took place, 
changes: R -v^ L. Simulations show that, whatever the initial configuration of 
scatterers, the particle always goes into a propagation phase in the direction 
of one of the lattice axes 0. An example is given in Fig.9. Here we show 
how propagation can be proved mathematically, and we derive the statistical 
properties of the dynamics of the particle. 

4.1 Definitions 

(i) The propagation of the particle on the triangular lattice always takes 
place in a strip which is defined a follows: a strip is a region of the triangular 
lattice bounded by two adjacent (parallel) lines, both oriented along one of 
the lattice axes. 

(ii) In accordance with the 1-D case, we define particle propagation as follows: 
a particle propagates in one direction in a strip if its motion is confined to the 
strip (defined by its initial velocity and its velocity after the first scattering 
event) where the particle visits any site not more than a fixed number of 
times. 

(iii) The particle's velocity vector at time t is denoted by C(t). 

4.2 The propagation theorem 

Theorem 2 : For any initial distribution of scatterers on a triangular lattice, 
a moving particle propagates in one particular direction through a strip on 
the lattice; this strip and the direction of propagation along it depend on the 
initial configuration of the scatterers at the origin (the initial position of the 
particle) and at three neighboring sites; these four sites form a parallelogram 
whose orientation determines the propagation strip. 

Proof : First we show that for each possible initial scatterers configuration, 
a propagation of the particle will occur in one of four strips, as sketched 
in Fig. 10 (propagation directions and strips are indicated by F(forwards), 
U(upwards), Di and D2 (downwards)). The strip in which the particle will 

^Although the spin notation is applicable to the rotators (with R =1 and L =|), the 
right (R) and left (L) characterization of the scatterers, as used in |^], is more convenient 
for our discussion of the triangular lattice. 
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propagate (F, U, Di or D2) depends on the initial velocity direction C{to), 
and on the initial state (R or L) of the four sites forming a parallelogram 
which contains C{to). There are four such parallelograms, hence four possible 
directions (see Fig. 10): 

®(2)(3)(6) ^ F, ®(2)@(7) ^ U, ®(2)(3)(8) D^, ®(2)@(5) D,. 
Then wc demonstrate that once the particle is inside a strip (at one of the 
sites indicated by (*) in Fig. 10), it never leaves this strip. Like in 1-D, the 
key ingredient of the ultimate unidirectional propagation of a particle on 
a triangular lattice fully occupied with flipping (R and L) scatterers, is the 
existence of a blocking mechanism. On the triangular lattice, this mechanism 
produces a blocking pattern based on a zig-zag path of length four (in lattice 
unit lengths), which contains parallel velocity vectors at even or at odd times 
(see examples in Fig. 11), and which prohibits particle motion more than one 
lattice unit length in the direction opposite to the propagation direction; as 
a consequence any site visited by the particle cannot be visited more than 
six times (see Corollary ^ below). The result is that the particle propagates 
in a strip in the direction naturally defined by the blocking pattern. 

Assume that a blocking pattern has been formed at the sites visited by 
the particle at times t,t + l,t + 2 and t + 3, so that the vectors C{t-\-) and 
C((t + 2)+), and C((t + 1)+) and C((t + 3)+) are parallel, respectively 
(see Fig. 11a). Then such a zig-zag trajectory can continue in two ways: 
either the particle continues directly its zig-zag motion at the next time step 
(see Fig. 12a), or it turns back (Fig.l2b). Obviously, the direct continuation 
of the zig-zag trajectory leads to propagative motion. So in order to prove 
propagation, we only have to consider the result of the turn back motion, and 
show that the particle nevertheless continues to move along the strip defined 
by the zig-zag direction along which it was moving from time t to time t-\-3. 
Now when the particle turns back, its trajectory creates a triangle (Fig. 12b) 
where from a zig-zag path emerges (Fig. 12c). Propagation then occurs as 
a consequence of the fact that backward motion at time t + 4, necessarily 
produces a blocking pattern, which forbids the particle to move more than 
one lattice unit length in the backward direction along the edge of the strip; 
as a result the particle will ultimately move one step further along the same 
strip in the zig-zag direction in which it was moving forward at time t + 3. 
Therefore any trajectory eventually forms a zig-zag of length 4, leading to 
propagation in the direction of the zig-zag path. One can verify that after the 
turning back of the particle, the trajectory forms a parallelogram (Fig. 12c). 
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Now the path of the particle through this parallelogram ends with three steps 
with velocity vectors: C((t + 7)+), C((t + 8)+), C((t + 9)+) that together 
form a zig-zag pattern of length 3 (see caption of Fig. 12). But these three 
velocity vectors form together with C((t+10)+) a new zig-zag of length 4, i.e. 
a new blocking pattern (such as in Fig. 12a), shifted one lattice unit length 
in the direction of propagation, with respect to the initial blocking pattern. 
The cases shown in the figures used for the demonstration are typical and 
all other cases are easily inferred from these typical situations, as the reader 
can verify. Consequently, we have shown that a blocking pattern always gets 
shifted by one lattice unit length along an edge of the strip in the direction 
of propagation, either after two time steps if the previous zig-zag motion is 
continued (Fig. 11a), or after seven time steps if the trajectory is turned back 
and a parallelogram of motion is formed (Fig. 12c). 

Thus, the particle will always propagate along a strip in the direction 
determined by the blocking mechanism in this strip. For an (arbitrary) given 
initial velocity of the particle, propagation can only take place in one of four 
strips: two with angles ±7r/3 (strips F and U) and two with angles ±7r (strips 
Di and D2) with respect to the initial velocity direction. The ratios of the 
occurrence of propagation along the four strips can be evaluated from Fig. 10 
by examining all possible paths leading to a strip (see examples in Fig. 11) and 
by giving weights q and (1 — g), to an R and an L scatterer, respectively. For 
q = \-i these ratios are 3:3:1:1 for the directions F, U, Di and D2 respectively. 

Corollaries follow from Theorem 2. 

Corollary 1 : The shortest time to create a blocking pattern is four time 
steps, and the longest time is ten time steps. 

Proof: The proof follows straightforwardly by inspection of Figs. 12a and 12c, 
where the zig-zag pattern of four velocities arranged in a W shape forms after 
4 and 10 successive displacements of the particle, respectively (as verified by 
counting the number of arrows after the initial state). 

Corollary 2 : The maximum number of visits to any site visited by the par- 
ticle on the triangular lattice is six, and the maximum number of passages 
along any link connecting two neighboring sites is five. 

Proof : Fig. 12c. shows the longest trajectory executed by the particle to 
move one unit length along the lattice in the propagation strip. When the 
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scatterer configuration is such that the particle is again forced into a turn 
back motion (as in Fig. 12b), when visiting the next site on tlie strip, the 
parallelogram of Fig. 12c repeats itself (shifted upside down by one lattice 
unit length in the propagation direction). Then the middle site of the par- 
allelogram, which has been visited three times (site (6) in Fig. 12c) is visited 
three more times, yielding a total number of six visits. The number of visits 
to a site is then easily evaluated by counting the number of incoming arrows 
pointing to the site. For instance, in Fig. 12c, there are three incoming arrows 
to site (6) , and. when at the next time step, the particle undergoes back- 
ward motion thereby producing again the parallelogram trajectory (shifted 
upside down to the right), the three arrows pointing towards site (3) are now 
pointing to site (6) , which yields a total number of six visits to site (6) . This 
is the maximum number of visits since the blocking mechanism prevents the 
particle from ever coming back to this site. Similarly the number of passages 
along a link between two neighboring sites is evaluated by counting the num- 
ber of arrows on that link. In the above example, which shows the typical 
case where the particle undergoes the longest possible trajectory to proceed 
in the propagation strip, the number of passages along link (3)-(6) is 3-1-2 = 5. 

Corollary 3 : All visited sites located on one edge of the propagation strip 
are in the same state (R or L), and all sites located on the other edge are in 
the opposite state (L or R, respectively). 

Proof : The blocking mechanism produces a zig-zag pattern with alternating 
velocity vectors (see Figs. 11 and 12); as a result the states of the visited sites 
must be alternately R and L, or L and R. Consequently, the propagation 
dynamics triggers a reorganization of the states of the scatterers of the visited 
sites. 

4.3 Statistical properties 

We characterize the propagation process by the increase in the value of the 
coordinate of the particle along one edge of the propagation strip, as illus- 
trated in Fig. 13. We consider the situation where the particle visits site r for 
the first time, and we analyze how the particle proceeds from site r to the 
(next) site with coordinate r + 1 on the same side of the strip. Obviously, 
before arriving at site r + 1, the particle must go through a site with coor- 
dinate r' = r + ^ on the other side of the strip (see Fig. 13). At both site r 
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and site r + ^, one of two possiblities occurs: either the particle trajectory 
continues a forward zig-zag path or it turns back. Thus, the transition from 
r to r + 1 can take place in four possibile ways: 

(//) the particle continues a forward zig-zag path at both sites r and r + ^; 
(fb) the particle continues a zig-zag path at r and turns back at r + |; 
(bf) the particle turns back at r, and goes into a forward zig-zag path at 
r + -■ 

{bb) the particle turns back both at r and at r + |. 

So, in case (//) the transition from r to r + 1 occurs in two forward time 
steps (see Fig. 12a), and obviously takes longer in the other cases. Indeed, 
in section ^]2|, we showed that if the particle turns back (see Fig. 12b), it 
takes seven (6+1) time steps before it performs a zig-zag path again along 
the propagation strip (see Fig. 12c). Therefore, the transition takes eight 
times steps in case (/6) - (1+7) corresponding to one step forward and one 
back turn - and (7+1) time steps in case (bf), while fourteen time steps are 
necessary in case (bb) - (7+7) corresponding to two back turns. It is easy to 
verify that the occurence probabilities of the various cases are: q{l — q) for 
case (//), for case (/&), (1 — g)^ for case (&/), and (1 — q)q for case (bb). 

From the above results, we can write the equation for the first visit prob- 
ability at site r + 1 at time t as 

Pi(r + l,t) = [q^ + (l-qf]P,(r,t-8) 

+ g(l-g)[Pi(r,t-2) + Pi(r,t-14)]. (4.1) 

It follows from this expression that the average time it takes the particle to 
propagate one lattice unit length along a side of a strip is (t) = 8 [g^ + (1 — 
g)^] + 16g(l — q) = 8 time steps, so that the average propagation velocity 
reads 

(c) = A = - , (4-2) 

(in lattice unit lengths per time step), independently of q and of the direction 
of propagation. 

When the particle is propagating in a strip, it takes an even number of 
time-steps to move from one site to the next on the edge of the strip in the 
forward direction (e.g. from r to r + 1 via r' in Fig. 13) because of the blocking 
zig-zag pattern.^ Therefore when the particle arrives for the first time at 

^"^A backward displacement (in the direction opposite to the propagation direction) can 
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site r + 1 (see Fig. 13), site r cannot have been visited more than four times 
(as can be checked with Fig. 12); so we can also write Eq. (^4.1|) as 



Pi(r + 1, t) = g(l - g)Pi(r, t - 2) + ^ P„(r, t - 2) 



(4.3) 



a=2 



where, as can be verified using Fig. 13 to explore all possible paths. 



P2(r,t-2) 
P3(r,t-2) 
P4(r,t-2) 



(l-g)2Pi(r,t-8), 

g2Pi(r,t-8), 

(l-g)gPi(r,t-14) 



(4.4) 



Now we define the total probability P(r, t) = J2a=i Pair, t), which, with 
is given by 



P(r,t) = P,(r,t) + [q^ + {l-q)'']Pi{r,t- 6) + q{l-q)Pi{r,t- 12). (4.5) 



In (|0| ), we use equation (^TT]) for Pi(r, t') with t' = t,t — Q,t — 12, and in 
the result, we combine the various terms to obtain a closed equation for the 
total probability 

P(r,t) = g(l - g)[P(r - l,t - 2) + P(r-l,t- 14)] (4.6) 

+ [q^ + {l-qf]P{r-l,t-8). (4.6) 

We observe that P{r,t) obeys the same equation as Pi(r, t), Eq.( ^yTD , and it 
can be shown ||^ that its continuous limit has the same structure as Eq. ( |3.43| ). 



4.4 Long-time behavior 

To investigate the long-time behavior of the particle motion, we proceed 
along the same lines as in the one-dimensional case. We start with the 
Laplace- Fourier transformation of Eq.( [4.1|) , which yields 



m=l 



only take place along the edge of the strip and takes one time step since, because of the 
blocking mechanism, it cannot exceed one lattice unit length. 
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g2 + (1 - qf 



e'^fkis) - E = 6 - 2m) 



g(l -g) 



m=l 
6 



+ 



(1 + e^2«)/,(s) - e^^shit = 12 - 2m) 



m=l 



(4.7) 



where fkif) denotes the space- Fourier transform of Pi(r, t), and /fc(s) its 
Laplace transform. In the hmit = 0, we obtain 



/fc=o(-s) 



9{s) 



(4.8) 



with 



5]e2-7,=o(14-2m) 

.m=l 



- g(l - q) 



Ee""7fc=o(12-2m) 

.m=l 

e2"^7fc=o(6 - 2m) 



.m=l 



and 



^?(sj 



,14s 



g(l-g)e^^^-[g^ + (l-g)2]e^^-g(l-g). 



(4.9) 



(4.10) 



To solve Eg. (|4.8| ), we need to know the "initial conditions", i.e. the values of 
fk(t) at t = 2m, for m = 0, ■ ■ ■ , 6. The values of /(r, t) at t = 0, ■ ■ ■ , 12 are 
obtained from the particle dynamics discussed in the previous section and 
are given in Appendix B, along with their corresponding Fourier transforms. 
We insert these values in ( [4.9| ), and, since we are interested in the long-time 
behavior, we consider the limit fk=o{s — » 0) by expanding fk=o{s) to lowest 
significant order in s; the result is 



lim/,=o(^) = l 



which, by inverse transformation, yields 

;im/,=o(t) = ^ 



(c). 



(4.11) 



(4.12) 



Thus, as for the one- dimensional case, we find that in the long-time limit, 
fk=o(t — >• oo) = limt^oo -Pi('^! ^) is equal to the average propagation ve- 
locity, with the difference that, in the triangular lattice, (c) is independent 
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of q. Furthermore it can be shown |^ that the long-time solution of Eq. ( |4.1| ) 
is given by 



(7r)^/' 



exp 



2'yr 



(4.13) 



with (t) = r/{c) = 8r and 7 = 72pq. In Fig. 14 we present numerical simu- 
lations and show that our theoretical result (|4.13| ) is in agreement with the 
simulation data. 



5 Concluding comments 

1. By replacing spins j and | by I's and O's respectively, the initial state of 
the 1-D spin lattice can be read as the input string of the tape of a Turing 
machine whose control device is the particle; the controller performs oper- 
ations according to the rules of the flipping spins and writes the output on 
the tape as the final spin configuration also converted into I's and O's. Here 
the computation algorithm performs the binary addition of the initial string 
with itself followed by a "XOR" with 1 (which is equivalent to 1 plus 
a shift). In the triangular lattice, one can similarly interpret the state of 
the scatterers along the edges of the propagation strip as R/L ^0/1; then 
the algorithm performs the logical operations "AND" and "NAND" with 
alternately on each side of the strip, which transforms a random sequence of 
O's and I's into a periodic string ...01010101... . 

2. After the passage of the particle through the one-dimensional lattice, 
all the spins are in the state opposite to their initial state (t-vv-J.) with one 
lattice site shift in the direction opposite to the propagation direction. Con- 
sider that the the system is made periodic (i.e. on a circle) by identifying the 
first and last sites of the chain, and that the particle propagates clockwise 
on the circle; then there is a counterclockwise drift (two lattice sites back) 
of the initial spin configuration every two cycles, and the average propaga- 
tion velocity is given by c(g) = ^[332^ + 3_2(i_g) ]- Similarly, when, in the 
triangular lattice, one attaches the two beginning and two end sites of a strip 
crosswise, the particle undergoes periodic motion on the strip. 
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3. Propagation has also been observed in the square lattice fully occupied 
with flipping scatterers when the scatterers are distributed periodically over 
the lattice 0. However, not in all such cases does propagation occur during 
the time span of the simulations (< 10'' time steps). It remains therefore an 
open question whether propagation on the square lattice, periodically cov- 
ered with flipping scatterers, always occurs or only for a certain class of pe- 
riodic scatterer patterns. When it occurs, propagation seems always to take 
place in a strip (like on the triangular lattice), which shows similarity with 
the glider behavior in the Game of Life [0 . Although the precise blocking 
mechanism is not known for the square lattice, "immediate" propagation can 
occur, which would imply then the possibility of an "instantaneous" blocking 
on the square lattice. 

4. Propagation on the one- dimensional and triangular lattices occurs for 
scattering angles of ivr and ±27r/3, respectively. For those cases on the 
periodically occupied square lattice where propagation has been found, the 
scattering angle equals ±7r/2, while for the randomly occupied square lattice, 
the maximum angle scattering rule yields a(n) (average) value larger than 
71/2. On the other hand, a scattering angle of ±7r/3 on the triangular lat- 
tice leads to the motion of the particle on a honeycomb (hexagonal) lattice 
where no propagation seems to occur for flipping scatterers 0. One might 
therefore conjecture that propagation only occurs in lattices fully occupied 
with flipping scatterers, when the scattering angle is > 7r/2 in absolute value. 

5. Notice that in the development performed for the triangular lattice, 
we have not used the regularity property of the Bravais lattice. So the results 
should not be dependent of this property and should also be valid for the 
random triangular lattice such as the Delaunay lattice (see Fig.l), provided 
the scattering rule is set such that the particle is deflected with the largest 
possible angle (> 7r/2). However, because of the topological randomness of 
the Delaunay lattice, one cannot predict whether the trajectory of the par- 
ticle will remain unbounded. 

6. To the best of our knowledge, it is an open question whether particle 
dynamics in three dimensional lattices with flipping scatterers can exhibit 
propagation or other peculiar behaviors. 
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Appendix A 



Here we derive the expression for the propagation velocity of a particle 
moving in a 1-D lattice with flipping spins, directly from the particle dynam- 
ics described in section The average asymptotic propagation velocity is 
given by 



W=,|mCW, (A.l) 



with 



C(t)^m^, (A.2) 

and where At is the unit time step (At = 1, i.e. the limit in (A.l) will be 
taken as t —>■ oo). It follows by iteration from Eq.( p.l| ), that 

R{t) - i?(0) = Cit - 1) + C(t - 2)) + ■ ■ ■ + C{0) = J2 C(r) . (A.3) 

T = 

When we consider the velocity of the particle at two successive time-steps, 
five cases must be examined as described in Table 1 (see also Fig. 15). If the 
particle arrives at a site at time r — 1 with velocity C(r — 1) and leaves the site 
with velocity C(r), we observe that when taking the sum over the velocities, 
all C{t — l)'s and C(r)'s which do not cancel are those which correspond 
to particle displacements leading to a site visited by the particle for the 
first time. So the sum Y.tC{'t) is equal to the number Nf of displacements 
which, during time t, lead to a first visit. Now, since the particle undergoes a 
displacement at each time-step, the total number of displacements Nt during 
time t, is equal to t. Consequently, from (A.1)-(A.3), we have 

(c) = hm \ [R{t) - i?(0)] = hm \ 1] C(r) = lim ^ ; (A.4) 

t^oo t t^oo t t^oo lytyt) 

and, since Nf{t)/Nt(t) is just Z^r -Pi('^) it follows that 

(c) = hm ^Pi(r,t) = hm h=o{t) = -\- ^ c{q), (A.5) 
where we have used (3.16). 
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The particle 


with 


and leaves 


The next site 


Fig. 


arrives at site 


velocity 


site r with 


is visited 


15 


r for the 


C(r - 1) 


velocity C(t) 


for the 


(a) 


first time 


C = +l 


C = +l 


first time 


(b) 


first time 


C = +l 


C = -l 


second time 


(c) 


second time 


C = +l 


C = +l 


first time 


(d) 


second time 


C = -l 


C = +l 


second time 


(e) 


third time 


C = -1 


C = +l 


second time 



Appendix B 



The "initial conditions" for the computation of the long-time behavior of 
the function fk=Q{t) = X]rPi(r, t) in section 4.4, are given by the following 
expressions 



f{r,t 


= 0) 






f{r,t 


= 2) 




- q) Sr,l , 


f(r,t 


= 4) 




- qf Sr,2 , 


f{r,t 


= 6) 


= q'{l 


- qf Sr,3 , 


f{r,t 


= 8) 


= q\l 


-qfSrA + [g' + (i-g)Vr,i, 


fir,t = 


= 10) 


= q'{l 


-qf6r,5 + 2g(l-g)[g2 + (l-g)Vr,2, 




.12) 


= q'{l 


-qfSr,e + 3q'{l-qf[q' + {l-qf]5r,s 



and correspondingly, in Fourier transform, by 



^(0) 


= 1, 






fk{2) 


= 


- q) e'' , 




^(4) 


= q\l 


-qfe^'^ 


•) 


fk{Q) 


= q'il 






fk{S) 


= q\l 




+ \q' + {l-q)fe^\ 


fkilO) 


= q'{l 


- q)^ gSjfe 


+ 2q{l-q)[q' + {l-qf]e'^\ 


fk{12) 


= q%l 


- qf e^'^ 


+ ^q\l-qf[q^ + {l-qf]e'^' 
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gure captions 



Fig. 1. Propagation in a random Delaunay lattice, where the particle 
arriving at a site of the random lattice is deflected over the largest pos- 
sible angle there, either to the right or to the left, depending on the R 
or L nature of the scatterer. Arrows illustrate particle displacements. 
The shaded shows the propagation strip. 



Fig. 2. Scattering on the one-dimensional lattice illustrating the proof 
of Theorem 1 (see text); the blocking patterns are framed with dotted 
squares. 



Fig. 3. Example of spin reorganization in 1-D lattice. The broken line 
is the particle trajectory. Black dots are sites with spin down, and open 
squares (shown only for the initial and final configurations) are sites 
with spin up. The spin reorganization is well observed by comparing 
the top row with the bottom row: • □ (with one site shift to the left). 



Fig. 4. Propagation velocity in 1-D lattice: for g = 1 (all spins initially 
up) (c) = 1 (upper line): q = 0.5 (for two different spin configurations; 
middle lines) (c) = 1/2; and g = (all spins initially down) (c) = 1/3 
(lower line). 



Fig. 5. fk=o{t) as a function of t for g = 0.1. Simulation results (dots) 
and theory, Eq. (|3.13|) (lines). 



Fig. 6. Space-time evolution of Pi{r,t) for q = 0.5. The x-axis denotes 
space. 



Fig. 7. The probability Pi(r, t) as a function of r at t = 500. Simula- 
tion data (black dots), binomial expression ( |3.23| ) (open circles), and 
Eq.dpl) (curve). 
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Fig. 8. Situations on the 1-D lattice leading to second visits, (a) and 
(b), and to third visit, (c). 

Fig. 9. Typical example of propagation pattern in triangular lattice. 
The initial position of the particle is shown by a full diamond; the 
propagation direction is along the upper left axis of the lattice. The 
coordinates are marked in lattice unit lengths. 

Fig. 10. Propagation strips on the triangular lattice. The initial ve- 
locity of the particle, C(to), is shown as the heavy arrow. When the 
particle arrives at one of the sites marked (I) , it is trapped in one of 
the propagation strips bounded by parallel heavy solid lines. The four 
propagation directions are indicated by arrows and capital letters (F, 
U, Di, D2). 

Fig. 11. Examples of trajectories leading to propagation in one of the 
four directions (F, U, Di or D2). Light arrows indicate the successive 
velocity vectors of the particle after the initial state {C{to), heavy ar- 
row). 

Fig. 12. Formation of blocking pattern by direct forward zig-zag path 
(a), and by turn back motion (b) followed by zig-zag trajectory emerg- 
ing from parallelogram path (c). Light arrows indicate the successive 
velocity vectors of the particle. Vectors C{{t+7)+),C{{t+8)+),C{(t+ 
9)+), and C((t + 10)+) are along the links (2)-(l) , (t)-® , ©-(*) , 
and (*)-0 ' respectively. 

Fig. 13. Particle motion from site r to site r + 1 along the propagation 
strip via site r' = r + ^. The arrow pointing to the right shows the 
propagation direction. 

Fig. 14. Long-time behavior in triangular lattice. Pi(r, t) as a function 
of r after 7616 time-steps: numerical data from simulation measure- 
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ments (black dots) compared with analytical result ( |4.13| ) (full curve). 

Fig. 15. Illustration of the five cases discussed in Appendix A and de- 
scribed in Table 1. 
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